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ABSTRACT 

In binary radio pulsars with a main-sequence star companion, the spin-induced 
quadrupole moment of the companion gives rise to a precession of the binary orbit. 
As a first approximation one can model the secular evolution caused by this classical 
spin-orbit coupling by linear-in-time changes of the longitude of periastron and the 
projected semi-major axis of the pulsar orbit. This simple representation of the pre¬ 
cession of the orbit neglects two important aspects of the orbital dynamics of a binary 
pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, due 
to the anisotropic 1/r 3 nature of the quadrupole potential. Secondly, the long-term 
secular evolution of the binary orbit which leads to an evolution of the longitude of 
periastron and the projected semi-major axis which is non-linear in time. 

In this paper a simple timing formula for binary radio pulsars with a main- 
sequence star companion is presented which models the short-term secular and most 
of the short-term periodic effects caused by the classical spin-orbit coupling. I also give 
extensions of the timing formula which account for long-term secular changes in the 
binary pulsar motion. It is shown that the short-term periodic effects are important 
for the timing observations of the binary pulsar PSR B1259-63. The long-term secular 
effects are likely to become important in the next few years of timing observations of 
the binary pulsar PSR J0045-7319. They could help to restrict or even determine the 
moments of inertia of the companion star and thus probe its internal structure. 

Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045- 
7319 since the analysis given in the literature is based on an incorrect expression for 
the precession of the longitude of periastron. A lower limit of 20° for the inclination 
of the B star with respect to the orbital plane is derived. 


Key words: pulsar timing binary pulsars - classical spin-orbit coupling - pulsars: 
individual: PSR .10045-7319, PSR B1259-63 


1 INTRODUCTION 


tract the maximum possible information from pulsar tim¬ 
ing observations, one needs an appropriate model (timing 
formula) for transforming each measured topocentric TOA, 
tobs, to the corresponding time of emission, T, measured in 
the reference frame of the pulsar. Various timing formulae, 
particularly for relativistic binary pulsars, have been devel¬ 
oped to describe radio-pulsar timing observations. 


Timing observations of pulsars, i.e. the measurement of the 
time-of-arrival (TOA) of pulsar signals at a radio telescope, 
is one of the few high-precision experiments in astronomy 
and therefore has a wide-ranging field of interesting applica¬ 
tions (Bell 1996). The first evidence for the existence of grav¬ 
itational waves as predicted by Einstein’s theory of gravity 
(Taylor & Weisberg 1989) and the first discovery of extraso¬ 
lar planets (Wolszczan & Frail 1992) are just the most strik¬ 
ing examples for the achievements of high-precision pulsar¬ 
timing observations. Approximately 10% of the known pul¬ 
sars are members of binary star system, i.e. in orbit around 
a white dwarf, neutron star or main-sequence star compan¬ 
ion. Timing observations of these binary pulsars is a power¬ 
ful tool to study various physical and astrophysical effects 
related to binary star motion and stellar evolution. To ex- 


With the discovery of PSR B1259-63 during a high fre¬ 
quency survey of the Galactic plane by Johnston et al. (1992) 
the first radio pulsar with a massive, non-degenerate com¬ 
panion was found. PSR B1259-63 is a 48-ms pulsar in a 
highly eccentric orbit with the main-sequence Be star SS 
2883. The second known radio pulsar with a massive, non¬ 
degenerate companion is PSR J0045-7319, discovered in a 
systematic search of the Magellanic Clouds for radio pulsars 
(McConnell et al. 1991, Kaspi et al. 1994). Some of the pa¬ 
rameters of these two main-sequence star binary pulsars are 
listed in Table 1. For both binary systems significant devi¬ 
ations from a Keplerian orbit have been detected which are 
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B1259-63 

J0045-7319 

P b (days) 

1237 

51.17 

x (sec) 

1296 

174.3 

e 

0.870 

0.808 

to (deg) 

138.7 

115.3 

i (deg) 

36 or 144 

44 or 136 

m, (M 0 ) 

~ 10 

8.8 ± 1.8 

Rt ( Rq ) 

~ 6 

6.4 ±0.7 


tagnon 1990). Finally, A so describes the Shapiro delay in 
the gravitational held of the Sun (Shapiro 1964). 

For binary pulsars the timing formula (j2j) has to be ex¬ 
tended by terms representing orbital motion and light prop¬ 
agation effects in the binary system (Blandford & Teukolsky 
1976, Damour & Deruellc 1986, Damour & Taylor 1992): 

T = fobs + A c — D / f 2 ± AflQ ± Abq + Asq 

(3) 

+A r + A e + A s , 


Table 1. Parameters of the main-sequence star radio binary pul¬ 
sars (Johnston et al. 1994, Kaspi et al. 1994, Bell et al. 1995). Pb, 
x, e, to, i, m * and it, denote respectively the orbital period, the 
projected semi-major axis, the eccentricity of the orbit, the lon¬ 
gitude of periastron, the inclination of the orbit with respect to 
the line-of-sight, the mass of the companion main-sequence star 
and its radius. 


most easily explained by classical spin-orbit coupling (Lai et 
al. 1995, Kaspi et al. 1996, Wex et al. 1997). Due to their 
high proper rotation the main-sequence star companions of 
PSRs B1259-63 and J0045-7319 show rotational deforma¬ 
tion and thus give rise to a gravitational quadrupole held. 
As a result of this a coupling between the orbital angular 
momentum and the spin of the companion takes place and 
leads to a precession of the binary orbit. 

In this paper I will shown that the present timing for¬ 
mulae represent only crude approximations to the orbital 
dynamics caused by the classical spin-orbit coupling and 
that there is a need for a new timing formula to model the 
timing observations. Before I focus on the construction of a 
new timing formula for main-sequence star binary pulsars a 
short presentation of the most important timing models is 
given. 

In a simple spin-down law the pulsar proper time is 
related to the phase, M, of the pulsar by 

M{T) = No + vT + \vT 2 + , ( 1 ) 

where v, v, and u are the rotation frequency of the pul¬ 
sar, its hrst and second time derivative, respectively (spin 
parameters). 

For a single pulsar the timing formula includes terms 
related to the position (a, 5), proper motion (/r a , Us) and 
parallax (-7r) of the pulsar. Moreover, it contains terms re¬ 
lated to relativistic time dilation and light propagation ef¬ 
fects in the solar system and also corrects for propagation 
effects in the interstellar medium (Backer 1989, Taylor 1989, 
Doroshenko & Kopeikin 1990): 

T = tobs + A c — D/f 2 

( 2 ) 

+A fl0 (a,5, /j, a ,ps, tt) + Abq + A S0 (a,<5). 

Ac corrects for the offset between the observatory clock and 
the ‘Terrestrial Dynamical Time’ represented by the best 
terrestrial standard of time. D/f 2 corrects for the disper¬ 
sive delay in the interstellar medium at the (barycentric) 
frequency fb where D is proportional to the column den¬ 
sity of free electrons between the pulsar and the observer. 
A_rq describes the so called Roemer delay, a change in the 
time of flight of the radio signal caused by the motion of 
the observer in the solar system reference frame. Abq rep¬ 
resents the transformation between ‘Terrestrial Dynamical 
Time’ and ‘Barycentric Dynamical Time’ (Fairhead & Bre- 


where the major effect is the Roemer delay Ar which de¬ 
pends on the orbital motion of the pulsar and the orientation 
of the pulsar orbit with respect to the line of sight. If the 
binary motion is purely Keplerian then the Roemer delay 
depends on 5 (Keplerian) parameters: 

• Pb , the orbital period of the binary system, 

• x = a p sin i/c, the projected semi-major axis, 

• e, the eccentricity of the orbit, 

• to, the longitude of periastron, 

• Tq, the time of periastron passage. 


a v is the semi-major axis of the pulsar orbit, i is the incli¬ 
nation of the orbital plane with respect to the line of sight, 
where i = 90° implies edge on, and c is the speed of light. 
The Roemer delay caused by the Keplerian motion of a bi¬ 
nary system is given by 

A r = a;[(cos U — e) sin a; + (1 — e 2 ) 1 ^ 2 sin U cos to] , (4) 


where U, the eccentric anomaly, is related to time, T, by the 
Kepler equation 


U — e sin U = 2n 


T-T 0 

Pb 


(5) 


Soon after the discovery of PSR B1913+16 (Hulse & 
Taylor 1975) it was clear that a pure Keplerian timing 
model is not appropriate to analyse the timing observation of 
this 7.8-hour orbital-period binary pulsar. For this purpose 
Blandford & Teukolsky (1976) derived a timing model (BT 
model) which contains the largest short-period relativistic 
effect, the ‘Einstein delay’ Ab, a combination of special- 
relativistic time dilation and gravitational redshift. They 
also included secular drifts of the main orbital parameters 
by following linear-in-time expressions: 


Pb = PbO + Pb (T — To) , 

(6) 

x = xo + x(T - T 0 ), 

(7) 

e = eo + e(T — To), 

(8) 

io = uiq + co(T — To). 

0) 


Based on a remarkably simple analytic solution of the post- 
Newtonian two-body problem (Damour & Deruelle 1985) 
Damour & Deruelle (1986) derived an improved timing for¬ 
mula (DD model) for relativistic binary pulsars. The DD 
model differs from the BT model in two ways: it contains 
the Shapiro delay As, which is of particular importance for 
i close to 90°, and it allows for periodic effects in the or¬ 
bital motion, e.g. in the BT model only the secular drift of 
the longitude of periastron is taken into account (equation 
(j|)), whereas the DD model describes both the secular and 
quasi-periodic changes in ui according to 

to = too + kA e (U), with to = 2nk/Pb, (10) 
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where 


A e (U) = 2 arctan 


l + ey / 2 

1-eJ 



( 11 ) 


and U is a function of T given by the solution of the (gen¬ 
eralised) Kepler equation 


U — e sin U — 2n 


( T -T 0 
\ Pi 


) 


h 

2 


( T -Tp 
V P b 



( 12 ) 


Pb accounts for any secular change in the orbital period, 
like the damping caused by tidal dissipation or the emission 
of gravitational waves. The post-Keplerian Roemer delay is 
given by 

A r = *{[cos U — e(l + 5 r )l sin u>(U) 

(13) 

+ [1 — e 2 (l + Sf) 2 ] 1 ' 2 sin Ucosu>(U)} , 

where u>(U) is given by equation ©■ The post-Keplerian 
parameters 8 r and 5f represent periodic post-Newtonian 
changes in the orbital motion, i.e. periodic changes of or¬ 
der ( v/c ) 2 where v is a typical orbital velocity of the binary 
star system. 

Damour and Taylor (1992) defined an improved version 
of the BT model (BT+ model) where equation (9) is re¬ 
placed by equation ©• The advantage of the BT+ model 
is that it contains the same number of parameters as the BT 
model but is able to account for quasi-periodic changes in 
U). 

To date, more than 50 radio pulsars are known which 
are members of a binary star system. The vast majority of 
these binary pulsars have a compact degenerate companion 
which is either a helium white dwarf or a neutron star. Tim¬ 
ing observations for most of these binary pulsars can be fully 
explained by the timing models above. In a first approxima¬ 
tion these models can be used for timing observations of the 
two main-sequence star binary pulsars, PSRs B1259-63 and 
J0045-7319, by fitting for the five Keplerian parameters, and 
for u! and x. However, this approximation models only the 
short-term secular changes of the binary orbit correctly. 

In this paper an improved timing model for binary radio 
pulsars with main-sequence star companions is presented. 
The new timing formula accounts for the short-term secu¬ 
lar precessional effects, for most of the short-term periodic 
orbital effects and for the long-term secular effects which 
are caused by classical spin-orbit coupling. In Section 2 a 
detailed investigation of the orbital dynamics of binary sys¬ 
tems with classical spin-orbit coupling is given. First, a sim¬ 
ple analytic solution is presented for the case that the orbital 
motion takes place in the equatorial plane of the massive 
companion. This solution will be helpful when developing 
the new timing formula. Then, the orbital dynamics of the 
general case, i.e. arbitrary orientation of the orbit with re¬ 
spect to the companion, is studied. In Section 3 it is shown 
how the orbital dynamics influences the timing observation 
of the main-sequence star binary pulsars PSRs B1259-63 
and J0045-7319. As a consequence the new timing formula 
is developed. In Section 4 the long-term validity of the old 
and new timing formulae is studied. Extensions which are 
quadratic in time and take into account long-term preces¬ 
sional effects in the binary orbit are presented. It is shown 
that observations of such long-term secular effects have the 
potential to probe the internal structure of the companion. 


In Section 5 I reinvestigate the orbital precession of PSR 
J0045-7319 since results presented so far in the literature 
were based on an incorrect formula for the precession of 
the longitude of periastron. In Section 6 the conclusions are 
given. 


2 THE ORBITAL MOTION OF 

MAIN-SEQUENCE STAR BINARY PULSARS 


The quadrupole of a main-sequence star companion is given 
by the difference between the moments of inertia about the 
spin-axis, 1%, and an orthogonal axis, /(. It is proportional 
to the mass of the companion, m*, the (polar) radius of the 
rotating star, R*, squared, and to the spin squared (Cowling 
1938, Schwarzschild 1958): 

I* — I* = rn*R*f2* = m*q , (14) 

where k is the apsidal motion constant and 12* = 

12*/(Gm*/i?*) 1 ^ 2 is the dimensionless proper rotation of 
the companion. 12* is the angular velocity of the compan¬ 
ion’s proper rotation. A main-sequence star of IOMq has 
k ^ 0.03, R, ~ 6 Rq. If the star is rotating at 70% of its 
break-up velocity, which appears to be a typical value for 
Be stars (Porter 1996), one finds for its quadrupole moment 

q = §fcR*f2* <, 0.35 R% = 7.6 x 10“ 6 AU 2 . (15) 


The spin-induced quadrupole moment of the compan¬ 
ion leads to an additional 1/r 3 potential term in the gravi¬ 
tational interaction between the two components which im¬ 
plies an apsidal motion and, when the spin of the compan¬ 
ion is not aligned with the orbital angular momentum, to a 
precession of the orbital plane. The general expressions for 
the rates of apsidal motion and orbital precession can be 
found in Smarr & Blandford (1976) and Kopal (1978), (see 
Section 4 in this paper). The expressions in Smarr & Bland¬ 
ford (1976) and Kopal (1978) are derived by averaging the 
orbital dynamics over a full orbital period, in order to get 
the secular changes in the orbit. This procedure, by defini¬ 
tion, neglects all short-term periodic effects. But, as will be 
shown in this paper, for main-sequence star binary pulsars 
with a long orbital period and a high eccentricity, like PSRs 
B1259-63 and J0045-7319, these short-term periodic effects 
are important. 

For a study of the short-term periodic effects one needs 
the orbital motion in full detail. In the centre-of-mass system 
the (Newtonian) orbital dynamics of a binary pulsar with an 
oblate companion star is given by the Hamiltonian (Barker 
& O’Connell 1975) 



GMn 


(1 + 5(5 [1 - 3(6 ■■»)’]) 


(16) 


where the linear momentum p is related to the linear mo¬ 
menta of pulsar and companion by p = p p = — p*. r is a 
vector pointing from the companion to the pulsar, r = |r|, 
n = r/r, M = m p + m* is the total mass of the system, 
H = m p m,/M is the reduced mass and s is the unit vector 
in direction of the spin of the companion. 

Before studying the full dynamics, I will investigate the 
special case of motion in the equatorial plane. Although the 
two known main-sequence star binary radio pulsars do not 
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orbit their companion in the equatorial plane the compara¬ 
bly simple solution of this problem will be helpful in devel¬ 
oping a new timing formula in the next section. 


2.1 The equatorial motion 


For motion in the equatorial plane (s ■ n = 0) the Hamilto¬ 
nian ( 0 ) reduces to 



The invariance of this Hamiltonian under time translation 
and spatial rotations implies the conservation of the (re¬ 
duced) energy, E = Ti±/n, and the (reduced) total angular 
momentum, J = r x p//i. If one introduces polar coordi¬ 
nates, r = r(cos/, sin/, 0 ), and makes use of the conserved 
quantities one finds the following equations of motion: 


or. , 2 GM J 2 , GMq 
2.E -\ - — --— 

V* <V»0 


(18) 


/ = 


r 2 ’ 


(19) 


where J = |J|. These equations of motion can be solved 
to first order in 5 by a quasi-Keplerian trigonometric 
parametrisation (cf. Damour & Deruelle 1985). For bound 
orbits, E < 0, one finds 


2 n ~~Fr~ = U — et sin U , 

( 20 ) 

r = a(l — e r cos U ), 

( 21 ) 

/ — /o = (1 + k)A ef (U ), 

where A ef (U) is given by equation (pd|) and 

( 22 ) 

2t tGM 

b (-2E) 3 / 2 ’ 

(23) 

(1 , 2EJ 2 2 EqV /2 

6t ^ G 2 M 2 J 2 ) 

(24) 

GM / Eq\ 

a = -W ( 1+ 

(25) 

e r = (1 + S r )et = ^1 - et , 

(26) 

, 3 / GM \ 2 

k =2{—) q ’ 

(27) 

e/ = (1 + 6f)e t = (i - ^jf) et ■ 

(28) 

The motion of the pulsar is 


r P = -Jl r ( cos /> sin /> °) • 

(29) 


Therefore, if the orbital angular momentum and the spin of 
the companion are (nearly) aligned then the orbital motion 
of the binary system can be described to first order in q by 
a simple trigonometric parametrisation which is identical 
to the ‘quasi-Keplerian’ parametrisation given by Damour 
& Deruelle (1985) to solve the post-Newtonian two-body 
problem. Further, the timing observations of such a system 
can be fully explained using the DD timing model. 



X 


Figure 1 . Definition of the angles 9, tp, and ip with respect to the 
equatorial coordinate system. / is the true anomaly and A* is the 
angle between the line-of-sight and the spin of the companion. 


2.2 The general case 

The full dynamics given by the Hamiltonian ( |lGl ) is a well 
known problem in the theory of Earth satellite motion. Due 
to the anisotropic nature of the quadrupole potential one 
does not expect any simple analytic solution as in the previ¬ 
ous section. Various methods have been developed to solve 
this problem (see e.g. Hagihara 1970, Roy 1978). To first or¬ 
der in q the dynamics given by equation can be solved 
by the method of the variation of the elements (osculating 
orbits). The following six elements are used to represent the 
osculating orbit: 

• a, the semi-major axis of the (relative) orbit 

• e, the eccentricity of the orbit 

• 9, the inclination of the orbital plane 

• <p, the longitude of the ascending node 

• ip, the longitude of periastron 

• M, the mean anomaly 


The angles are defined with respect to the equatorial plane 
of the oblate companion (see Fig. 1). 

Lagrange’s planetary equations for this problem can be writ¬ 
ten in the form 


da 2 977 

dt na dM ’ 

de (1 — e 2 ) 1 / 2 / _ 2 . 1/2 dn _ 977 \ 

dt na 2 e f dM dip J ’ 

dd _ _1_ f 977 977 \ 

dt na 2 (l — e 2 ) 1 / 2 sin 9 y dip d(p J ’ 

d4 _1_ 977 

dt na 2 (l — e 2 ) 1 / 2 sin# d9 

dip _ 1 (l -e 2 dn 977 \ 

dt ~ na 2 (l — e 2 ) 1 / 2 V e de C ° & d9 ) ’ 


(30) 

(31) 

(32) 

(33) 

(34) 


dM _ 2 / 1 - e 2 dll 977 \ 

dt na 2 ( 2 e de + da J 


(35) 


The relevant parts of the disturbing function 1Z = 1Z S + TZp 
are 
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Ti s = 


Tip = 


qn 2 / 3 . 2 a \ /, 2 x- 3/2 

*r (0/r)> 


(1 — e 2 ) 3 / 2 


+ — sin 2 0 cos 2 (ip + f) 


(36) 

(37) 


TZs and Tip are first-order secular and short-period parts 
respectively of the disturbing function li. 

The secular perturbations of the first order are obtained 
by putting 7?. = 77. s in equations (bd) to (B3). The result is 


A a s = 0 , Ae s = 0 , A 9 S = 0 , 
A<p s = —Q cos 0 O n 0 t, 

Aip s = Q (2 — 3 sin 2 0o) not, 


A M s = 


■(> 


1 + Q (1 - ^ sin 2 0 O ) (1 - eo) 


1/2 


not, 


where the definition 


(38) 

(39) 

(40) 

(41) 

(42) 


was used. The zero subscript indicates evaluation at the ini¬ 
tial condition, no is related to the initial semi-major axis, 
a 0 , by 


n 0 a 0 — GAT . 


(43) 


To simplify the representation of the first-order short- 
periodic perturbations I define the functions 

3 • 2 , 

— sir 

2 


(-!-■•) 


1 - yj- J sin / + | sin 2/ + £ sin 3/ 


Ce,ej = (l - | sin 2 


(44) 


2 \ 2 
1 + ^- ) cos / + | COS 2/ + -^ cos 3/ 


and 

S^-j{a 0, «1, Q2, 03, 04, as) = 

5 

Qo sin( 2 i/> — /) + ^ on sin( 2 i/) + mf ), 

m=i ( 45 ) 

C^ f (ao, a i, Q 2 , Q3, 04 , as) = 

5 

ao cos( 2 tp — f) + ^ an sm(2tp + mf), 


and 

W e j = f — M + e sin / . 


(46) 


To derive the first-order short-period perturbations, the dis¬ 
turbing function 7 i in equations ( |3o| ) to (^i|) is replaced by 
Tip. Integration of the resulting equations leads to the follow¬ 
ing expressions for the six elements (Kozai 1959, Fitzpatrick 
1970): 

2eoao , 

A a p = Q- -j {6'eo.flo;/ 


1 -el 


+ sin 2 0 o 1 


( <=o 12 +3<=o 

2+3ej( 

12+3e= 

3en e oAl 

l 16’ 16 ’ 

4e 0 ’ 

16 ’ 

8 ’ 16 J j 


(47) 


(48) 


Ae p — Q {Ceo,0o;/ 

+ s *(& 4+lle 0 5e 0 28+17e 0 3e 0 eg\\ 

' Wo;/ l 16’ 16 ’ 4 ’ 48 ’ 8 ’ 16 J j 

A0 p = Qcose 0 sm0 0 C UJO .f(0,f,^f,0,0) , ( 49 ) 

A^ = Qcos0o{-IF/ + ^ o ; / (O,f,i,a,O,O)} , (50) 

AiA P = £ { ( 2 - ^ sin 2 0 O ) e 0 Wf + S eo ,g 0 , f 
-eo Sto,f (0,f,i,f,0,0) 

+ sin 2 0 ( 


(51) 


( 

4-15e§ 

5e 0 

28+19e^ 

3eo 


y 16 ’ 

16 ’ 

4 ’ 

48 ’ 

1 8 ’ 

16 Ji 


* U - n^~ e o ) 1/2 ro 

AA 4 P — —y-{a eo ,e o; / 

eo 

+ sin 2 0 O S’vo,/ 


196 —7ej^ 3ej 
336 ’ 8 


-*)}■ 


(52) 


The mean and true anomaly are connected by the equa- 

(53) 


tions 

A4 = U — e sin U , 

'1 + eX 1 / 2 


/ = 2 arctan 


(rH) 


tan - 


U 


(54) 


The distance between the main-sequence star and the pulsar 


r = a(l — ecos U) = 


a(l — e 2 ) 
1 + e cos / 


(55) 


The position vector of the pulsar at the time t with respect 
to the equatorial coordinate system is 


( cos ij> — sin <p 0 
sin tp cos cj> 0 

0 0 1 


cos (tp + f) \ 
sin(i)> + /) cos0 J . (56) 

sin(i)> +/) sin 0 / 


where the (variable) parameters a, e, 0, (j), ip, M are given by 

€(*) = Co + [A Ut) - A£,(0)] + [A4(f) - A&(0)] (57) 

(£ stands for one of these six parameters). To calculate A4 
one uses M = Mo + not, e = eo in equations ( |53| ) and (|54j) 
to obtain an approximated value for / which then is used in 
evaluating equations @ to (p^). 

As an example I substitute the orbital elements of the 
main-sequence star binary pulsar PSR B1259-63 (Table 1) 
into equations (p7| ) to (pi|), For the initial value of 0 (0o) I 
take 60°, which is a realistic value for this binary system (see 
Wex et al. 1997). The values for tp o and ipo can be estimated 
from lj, i, and A*: 


cos i — cos 0 cos A* 
sin 0 sin A* 
sin A* sin <j> 


sin(ui — ip) = — 
cos(u> — ip) = 


sin * 

sin 0 cos A» — cos 0 sin A* cos <p 
sin i 


(58) 

(59) 

(60) 


(Equation (p^) does not determine <p uniquely, since <p runs 
between 0 and 27r.) A* can be derived from optical observa¬ 
tions of the projected proper rotation of the Be star. John¬ 
ston et al. (1994) found A* ~ 30° (or 150°). For the ‘strength’ 
of the quadrupole I use q = 2 x 10 -6 AU 2 (fc ~ 0.01, cf. equa¬ 
tion (B). 
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The results for two full orbits of PSR B1259-63 are 
given in Fig. 2 and Fig. 3. Fig. 2 presents the changes of the 
four ‘periodic’ parameters, a, e,9, M, which take their initial 
value after every full period. Fig. 3 presents the changes of 
(f> and i!>. I call cj) and if> ‘secular’ parameters since they show 
both periodic and secular changes. It is obvious that for all 
parameters the major changes take place very close to the 
periastron passages, as expected from the 1/r 3 nature of the 
quadrupole potential. 


3 TIMING MODELS FOR MAIN-SEQUENCE 
STAR BINARY PULSARS I. SHORT-TERM 
PERIODIC EFFECTS 


Let Ko = (0, — sin A,, cos A*) be the unit vector which indi¬ 
cates the direction of sight (see Fig. 1). The Roemer delay 
measured by an observer on Earth is then given by 


A r = - Ko ■ r p , 
c 

where c is the speed of light and 
to* 

Fp = "M r 


(61) 


(62) 


is the position vector of the pulsar originating in the centre 
of mass of the binary system. Using equation ( |56| ) leads to 

Afl = -jj ~ {—sinA*[sin < ()COs('i/> + /) 

+ cos 4> sin(i p + /) cos 9] + cos A* sin (ip + /) sin 9 } . 


As mentioned in the introduction, the simplest timing 
model for binary pulsars is the BT model where changes in 
the longitude of periastron, t o, and changes in the projected 
semi-major axis of the pulsar orbit, x, are assumed to be 
linear in time. The discussion in the previous section clearly 
showed that the (osculating) parameters of a binary system 
with an oblate companion do not change linearly in time (cf. 
Fig. 2 and Fig. 3). Thus one does not expect that the appli¬ 
cation of the BT model leads to a perfect fit. Fig. 4 shows 
the difference between equations (^)>(@),@ and the actual 
Roemer delay expected for the binary pulsars PSRs B1259- 
63 and J0045-7319. (Using the DD model instead leads to 
similar results). The typical precision in the measurement 
of the arrival time of pulsar signals is of the order of 100 /ts 
for PSR B1259-63 (Johnston et al. 1994) and a few ms for 
PSR J0045-7319 (Kaspi et al. 1994). For both pulsars the 
deviations given in Fig. 4 are larger than the error in the 
TOAs. For PSR B1259-63 it is more than a factor of ten. 
Therefore the BT model is only a very crude approximation 
for these two binary pulsars and there is the need for an 
improved timing model, for PSR B1259-63 in particular. 

To extract reliable information from timing observa¬ 
tions of main-sequence star binary pulsars one could con¬ 
struct a timing model that contains the full orbital dynamics 
given in the previous section. In the (unlikely) case that the 
orbital motion takes place in the equatorial plane of the Be 
star one can use the DD timing model as shown by the so¬ 
lution in section 2.1. For the general case one has to use the 
solution of Section 2.2 which leads to a rather complicated 
timing formula with a comparably high number of param¬ 
eters to fit for where some of these parameters are only in¬ 
directly related with observable effects. Therefore, given the 
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Figure 2. Calculated changes of the ‘periodic’ parameters of PSR 
B1259—63. I used the parameters of Table 1 (?.q = 36°), A* = 30°, 
9 0 = 60°, (j >o = 30°, ipo = 114°, and q = 2 X 10“ 6 AU 2 . 


limited number of TOAs and the finite size of their mea¬ 
surement errors, one sees that in general a timing formula 
based on the equations of Section 2.2 is not a practical pro¬ 
cedure. What one is looking for is a simple timing model 
which is a very good approximation to reality. In the ideal 
case the number of parameters should be the same as in the 
BT model. 
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Figure 3. Calculated changes of the ‘secular’ parameters of PSR 
B1259—63. (Parameters as in Fig. 2.) 


The main deficit of the BT model is the use of equations 
(i~) and (Ji|) to describe the precession of the orbit. Even for a 
purely equatorial motion (see Section 2.1) equation (^j) has 
to be replaced by ([HJ according to the DD timing model. 
If the orbit is tilted with respect to the equatorial plane the 
precession of the orbit leads to a change in the inclination 
of the orbital plane with respect to the line of sight. Fig. 5 
gives this change of i for PSR B1259-63 as a function of the 
time, t, and as a function of the true anomaly, /. The change 
of i is neither linear in t nor linear in /. But the assumption 
of linearity in / is obviously much closer to reality than the 
assumption of linearity in t. The same is true for the change 
of x which is a function of i. 

Therefore I define a new timing model which I call BT++ 
model. Analogue to the construction of the BT+ model by 
Damour and Taylor (1992) the BT++ model is defined by 
replacing equations (Q) and (fi) in the BT model by 


x = x 0 + £A e (U ), 

and 

x = 2n£/P b , 

(64) 

LO = LOq kA e (JJ ), 

u> = 2nk/P b . 

(65) 


A comparison of the expected Roemer delay and the Roe- 
mer delay as used in the BT++ model is given in Fig. 6. For 
both main-sequence star binary pulsars, PSRs B1259-63 and 
J0045-7319, the BT++ model is off by clearly less than the 
typical error in the TOAs. Only very close to periastron the 
deviations for PSR B1259-63 show a sharp peak of 200 (is, a 
value which is slightly larger then the typical measurement 
error. On the other hand, so far there are no timing obser¬ 
vation of this pulsar close to periastron. The reason is the 
occultation by the circumstellar disk which lasts from ~20 



orbital phase (t/P^ 



orbital phase (t/PfJ 


Figure 4. Roemer delay as used in the BT timing model minus 
Roemer delay as expected for PSRs B1259-63 (upper figure) and 
J0045-7319 (lower figure). For B1259-63 I used the same parame¬ 
ters as in Fig. 2. For J0045-7319 I used the parameters of Table 1 
(to = 44°), A* = 15°, 0 O = 30°, (j )o = 156°, ^o = 106°, and 
q = 2 X 10“ 6 AU 2 . 


days before until ^20 days after periastron (Johnston et al. 

1996) . 

I conclude that one should certainly use the BT++ 
model instead of the BT model to fit the TOAs of PSR 
B1259-63. The BT++ model has the same number of pa¬ 
rameters as the BT model, but is able to account for the 
fact that changes of the binary parameters happen mainly 
close to periastron. The BT++ model was already applied 
successfully to fit the TOAs of PSR B1259-63 (Wex et al. 

1997) . I expect a slight improvement of the residuals for PSR 
J0045-7319 when one uses the BT+-1- model. 

As concluded in Section 2.1, the correct timing model 
for equatorial orbits is the DD model. The DD model takes 
into account all the periodic effects of the (equatorial) orbital 
motion. One can now try to construct an even better timing 
model for main-sequence star binary pulsars, say DD+, by 
replacing equation (7j) in the DD model by equation (|64|). 
The result is a timing model which combines the advan¬ 
tages of the BT+-1- model in describing the precession of 
the orbital plane and the DD model in describing periodic 
orbital effects. The representation of the Roemer delay in 
the DD+ model contains one more parameter than in the 
BT++ model. The DD+ model has the same number of pa¬ 
rameters as the DD model. From Fig. 7 one sees that the 
DD+ model represents a nearly perfect fit for most parts of 
the orbit and close to periastron it is an improvement by 
a factor of 2 with respect to the BT+-(- model. At present 
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Figure 5. Calculated change of the inclination i of the orbit of 
PSR B1259-63 as function of the time, (upper figure) and of 
the true anomaly, /, (lower figure). (Parameters as in Fig. 2) 




Figure 6. Roemer delay as used in the BT++ timing model mi¬ 
nus Roemer delay as expected for PSRs B1259-63 (upper figure) 
and J0045-7319 (lower figure). (Parameters as in Fig. 4.) 



Figure 7. Roemer delay as used in the DD+ timing model minus 
Roemer delay as expected for PSRs B1259-63. (Parameters as in 
Fig. 2.) 



Figure 8. Roemer delay as used in the step-function timing 
model minus Roemer delay as expected for PSRs J0045-7319. 
(Parameters as in Fig. 4.) 

the measurement precision for the TOAs for PSRs B1259- 
63 and J0045-7319 does not allow to fit for the (full) DD+ 
model. 

Finally, the upper figure of Fig. 5 indicates that a step- 
function model (a> and x change discontinuously at each pe- 
riastron) will also give a good approximation to reality. In 
fact, Fig. 8 shows that a step-function model is clearly bet¬ 
ter than the BT model (lower figure of Fig. 4). On the other 
hand, comparison between the lower figure of Fig. 6 and 
Fig. 8 implies that a step-function model is clearly worse 
than the BT++ model, in particular for observations close 
to periastron. 


4 TIMING MODELS FOR MAIN-SEQUENCE 
STAR BINARY PULSARS II. LONG-TERM 
SECULAR EFFECTS 

So far only the short-term effects in the orbital motion have 
been dealt with. I have shown the advantage of the BT+-f 
(and DD+) model in taking into account the short-term 
periodic effects of the orbital motion. In the BT, BT+, DD, 
BT++ and DD+ model the secular changes in u> and x are 
assumed to be linear in time. This approximation will hold 
as long as there are only small changes in u> and x. In this 
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Figure 9. Definition of various angles in the total-angular- 
momentum reference frame. The invariable (X-Y) plane is per¬ 
pendicular to the total angular momentum J = L + S. The line- 
of-sight (Ko) is in the Y-Z plane. J is a conserved quantity and, 
if averaged over one full period, jL| and |S| are conserved quanti¬ 
ties. Thus ij, 9j, and 9 are conserved. The angles 4> and \F change 
linearly with time. 


section I will focus on the long term precession of the binary 
orbit and its influence on pulsar timing and will investigate 
the limits of the present timing models. 

In the following discussion I neglect periodic effects and 
focus only on the secular changes in the orbit of the bi¬ 
nary system caused by the spin induced quadrupole of the 
main-sequence star companion. The solution presented in 
Section 2.2 does not give the long term behaviour correctly. 
It does not take into account the change in the orienta¬ 
tion of the main-sequence star due to spin-orbit coupling. 
The change of the orientation of the main-sequence star is 
of order q and thus appears in the equations of motion at 
order q 2 which was neglected in Section 2.2. On long time 
scales the orientation of the main-sequence star changes by 
a comparably large amount and therefore the contribution, 
although of order q 2 , becomes numerically significant. The 
solution in Section 2.2 is perfectly suited for a discussion 
of periodic effects during a few orbital turns, but for the 
study of the long-term behaviour one should focus on the 
conserved quantities, which are the total energy and the to¬ 
tal angular momentum. The total angular momentum, J, is 
the sum of the orbital angular momentum, L, and the spin 
of the main-sequence star, S. On average, over one full pe¬ 
riod, the length of S, |S|, and L, |L|, are conserved (Barker 
& O’Connel 1975). Fig. 9 illustrates the resulting orbital 
dynamics. 

Averaged over one orbital period one finds for the 
change of <f> and <F (Smarr & Blandford 1976, Kopal 1978) 


$ = -nQ 


sin 9 cos 9 
sin 9 j 


= const. 


and 

'F = nQ ^1 — ^ sin 2 9 + i sin 29 cot 9jJ = const ., 
where 


( 66 ) 


( 67 ) 


p 3(J 3 * - II)/m, kR 2 M 

^ 2a 2 (l-e 2 ) 2 a 2 (l — e 2 ) 2 


( 68 ) 


The bar on top of the quantities indicates that they are av¬ 
eraged over a full orbital period (c.f. Section 2.2). For sim¬ 
plicity I will skip the bar on top of the averaged quantities 
for the rest of this section. Equations (^6|) and ( |67| ) can be 
derived directly from equations (jl9|) and (f!o|). 

For the inclination of the orbit with respect to the line 
of sight, i, and the longitude of periastron, one finds the 
relations 


cos i = cos ij cos 9j — sin ij sin 9j cos $ , (69) 

sin a; = [(sin 9j cos ij + cos 9j sinij cos <f>) sin T 

sm* (qU) 

+ sinzj sin $ cos *F), 


cos uj = ——: [(sin 9j cos ij + cos 9j sin ij cos $) cos 'F 
sin i 

— sin ij sin $ sin <F) 


(71) 


(see Fig. 9 for the definition of ij, 9j, $ and T). The angles 
ij and 9j are conserved quantities. The angles <F and >F 
evolve linearly in time, t, 

$ = $o + $(4 — to) and *F = To + 4>{t — to) ■ (72) 

Equations ( |69| ) to ( fr2| ) give the full (secular) evolution of the 
projected semi-major axis, x = a p sin i/c, and the longitude 
of periastron, ui. This evolution is clearly non-linear in time 
since changes in >F couple in a complicated way with changes 
in $ to produce the secular changes in i and u). 

A first approximation of this non-linear behaviour can 
be given by ^ 


x(t) ~ xo + xo(t — to) + (f — to) 2 , 

(73) 

and 


u(t) ~ cuo + ch 0 (f - to) + §iio(t — to) 2 . 

(74) 


The quantities xo, xo, Coo, u>o, if measured in timing obser¬ 
vations, contain information about the orientation and the 
quadrupole moment of the companion star. 

In general the orbital angular momentum is much larger 
than the spin of the companion star and thus < 1. In 
case of the binary pulsars PSRs B1259-63 and J0045-7319 
|S|/L ~ 0.1 and 9j ^ 6°. If 9j i,7r — * one finds after an 
(Laurent) expansion of equations (p9|) to ( fn| ) with respect 
to the small angle 9j: 

i = ij + 9j cos <F + 0(9 2 ), (75) 

cu = IF + <f> — cot i sin $ + 0(9 2 ), (76) 

and 


xo = nQx o cot i sin 9 cos 9 sin $0 + 0{9j ), (77) 

+ O(9 0 j ), (78) 


2 „2 . ( sin 9 cos 9 cos $0 

xo = —n Q xo cot 1 


u>o = nQ (1 — ^ sin 2 


sin 9j 

9 + cot i sin 9 cos 9 cos $ 0 ^ +0{9j) ,(79) 


T To include short-term periodic effects, as discussed in the previ¬ 
ous section, one has to replace xnjt — to) by £,A e (U) and iho(t —to) 
by kA e (U ); see equations (f i4 [), (Ibol). 
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ido = n 2 Q 2 cot i 


sin 2 9 cos 2 9 sin $0 
sin 9j 


+ O(9 0 j) . 


(80) 


While doing these expansions it is important to keep in mind 
that and the leading term of are of order 9J 1 but < f? + 'F 
is only of order 9 j. A fact which has been overlooked by 
Smarr and Blandford (1976) and Lai et al. (1995) leading to 
a wrong result for loo- 

Since Q is a quantity which is not well known, it is 
not possible to extract the angles 9 and $0 just from the 
measurement of io and u>o- But, if the masses, and therefore 
i, are known within a certain accuracy (which I assume for 
the following considerations) the measurement of io and ujo 
restricts possible values of 9 and $0 to a small region in 
the 9-&o plane by dividing equation ( friil ) by equation ([n]), 
which leads to 


UJOXO . , , 

-sm <t >0 — cos <t >0 = 

io 


1 + 3 cos 29 

2 cot i sin 29 


(81) 


The sign of io or u>o will help to exclude further regions in 
the 6 > -<f»n plane. Since the right-hand sides of equations ( ]T7[ ) 
and (179) have to be positive or negative depending on the 
sign of io and u>o respectively (c.f. Fig. 11). 

If one is able to measure either x or id, in addition to i 
and lj, one can make use of the following restrictions on $0 
and 9j: 

io tan i . . , , „ , on . 

- = — sm i&o tan $0 9j (82) 

*0*0 

or 


Xq tan i 
idQX% 


= sin $0 9j . 


(83) 


If one is able to measure i, tjj, x, and Co for a binary pulsar 
with a main-sequence star companion, then 


id 0X0 
X 


— tan $0 


(84) 


will give + 0 . Equation (^2|) or (js^) will then give 9j and 
equation (| 8 lj) will give 9. Therefore all angles in the geome¬ 
try of the binary system are determined and one can further 
get Q from equation (|77|) or (Jt 9[) and the spin of the compan¬ 
ion using |L|, 9j and 9 (see Fig. 9). This allows to determine 
the moments of inertia J 3 and I{ of the companion, which 
are related with its internal structure. 

For PSR B1259-63 the changes of i and c d are typically 
of the order of one second of arc per orbit (see Fig. 5). Since 
the orbital period is 3.4 years, the quadratic-in-time changes 
of these angles will be absolutely negligible for the next few 
decades. 

For PSR J0045-7319 the situation is different. The 
changes in x and u) are two orders of magnitude larger than 
for PSR B1259-63. Fig. 10 shows the result of fitting for 
simulated pulse arrival times for PSR J0045-7319. The first 
figure (10a) gives the pre-fit residuals for a BT model that 
leads to a good fit for the first few orbits. After four years 
the model is off by about 40 ms. Fitting for the whole time 
span of 1500 days using the BT model one finds the post-fit 
residuals given in the second figure (10b). Most of the de¬ 
viations of Fig. 10a are absorbed in the spin parameters by 
changing them according to 


A P = -2.1 x 1CT 10 s, A P= -3.3 x 1CT 18 . (85) 


The residuals are in the order of the present measurement 



t (days) 



t (days) 



Figure 10. Simulated pulse arrival times for PSR J0045-7319 
over a time span of 30 orbits (~ 4.2 years). Top: pre-fit residuals 
for a BT model which leads to a good fit for the first few orbits. 
Middle: post-fit residuals after fitting for spin and orbital param¬ 
eters, including x and lj. Bottom: post-fit residuals after fitting 
for spin and orbital parameters, including x and lj, and second 
derivatives, x and lj. I used the values of Table 1 (i = 44°), 
6 = 30°, 4> 0 = 80°, k = 0.01, and ft* = 0.5. |S|/|L| « 0.1 and 
thus 0j « 3°. 


precision. For future observations the BT model (and there¬ 
fore the BT++ and DD+ models) will fail to explain the 
observations and one has to fit for higher derivatives in x 
and lj. The result of such a fit is presented in the last figure 
( 10 c). 

Finally, it should be mentioned that fitting for a 7+ 
instead of x and id improves the residuals only marginally 
and gives a P;, which is two orders of magnitudes smaller 
than the one observed in this system. Therefore the non¬ 
linear drifts of x and u> cannot explain the observed Pb . 
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5 PSR J0045—7319 — EVIDENCE FOR A 
NEUTRON-STAR BIRTH KICK? 


The orbital precession of the binary pulsar PSR J0045-7319 
was seen as a direct evidence that the neutron star of this 
system received a kick of at least 100 km/s at the moment 
of birth (Kaspi et al. 1996). Since the theoretical analysis in 
this paper is based on the calculations of Smarr and Bland- 
ford (1976) and Lai et al. (1995) an incorrect formula for 
u), the precession of the longitude of periastron, was used 
(equation (1) in Kaspi et al. 1996). Therefore their limits on 
the angle between the spin axis of the B star and the orbital 
angular momentum, 25° < 6 < 41°, are incorrect and one 
has to reinvestigate the question whether this binary star 
system provides an evidence for a neutron-star birth kick. 

Using the values of Table 1 in Kaspi et al. (1996) one 
finds for PSR J0047-7319 


± 0*0 
X 0 


1.80 ±0.05. 


( 86 ) 


Co was corrected for the general relativistic contribution of 
0.004°/yr. The uncertainty in the masses leads to 


i = 44° ±5° (or * = 136° ±5°). (87) 

Equation (^I|) restricts the values of <f>o and 9 as shown in 
Fig. 11 and one finds as a lower limit on the inclination of 
the B star with respect to the orbital plane 


9 > 20° . 


( 88 ) 


This limit is only slightly smaller than the one given in Kaspi 
et al. (1996) and so does not change their major conclusion, 
i.e. that the binary system PSR J0045-7319 provides direct 
evidence for a neutron-star birth kick of at least 100 km/s. 
However, their conclusion that ±o > 0 implies 9 < 55° is 
incorrect. In principle one can have 6 up to 70°, although 
9 being close to 70° requires an unphysically high apsidal 
motion constant k for the B star. Numerical simulations 
show that 6 > 65° is excluded by the fact that present timing 
observations are still in agreement with a simple x-Co model. 
Similar arguments apply for the case 9 > 90°. 


6 CONCLUSIONS 

In this paper I have presented a timing formula for main- 
sequence star binary pulsars which takes into account most 
of the periodic variations along the orbit caused by the 
anisotropic nature of the 1/r 3 potential of the spin-induced 
quadrupole moment of the companion star. The new tim¬ 
ing formula contains the same number of parameters as the 
Blandford-Teukolsky timing formula. I have shown by nu¬ 
merical simulations, that the new timing formula leads to 
much better results in case of the long-orbital period binary 
pulsar PSR B1259-63 then the Blandford-Teukolsky or the 
Damour-Deruelle timing formula. Only very close to perias¬ 
tron the new timing formula shows deviations are slightly 
greater than the typical measurement error in the time-of- 
arrival of the pulsar signals. But so far there are no timing 
observations of PSR B1259-63 close to periastron, due to 
the eclipse of the pulsar caused by the circumstellar mate¬ 
rial. For PSR J0045-7319 these periodic variations are of 
the order of the measurement precision. 

I have given quadratic-in-time extensions of the timing 



Figure 11. Observational constraints on the geometry of the PSR 
J0045-7319 binary system. I used i = 44° ± 5° (for i = 136° ± 5° 
one has to replace 4>o by +q — 180°). The grey areas are excluded 
because x > 0. The region between the two dashed curves is 
excluded because ojq > 0. (To plot the dashed curves I used i = 
44°). The solid curves border the narrow region of possible values 
of 8 and ( l'o. The width of these regions is dominated by the error 
in i. The regions within the dotted curves are ruled out since 
the proper rotation of the B star (projected value: 113+10 km/s) 
does not exceed the break-up velocity (~420 km/s), see Lai et al. 
1995. These regions do not lead to any further constraints on 8 
and 4>. 


formula which account for long-term secular changes in the 
orientation of the binary-pulsar orbit. In particular for the 
binary pulsar PSR J0045-7319 these extensions might be 
important in the next few years of timing observations de¬ 
pending on the, so far unknown, orientation of the B star 
spin and the total angular momentum of the binary system. 
I have concluded that the measurement of these long-term 
secular effects has the potential to probe the internal struc¬ 
ture of the companion. 

Finally I have reinvestigated the classical spin-orbit pre¬ 
cession of the binary pulsar PSR J0045-7319 since the the¬ 
oretical analysis of this binary system given in Lai et al. 
(1995) and Kaspi et al. (1996) is based on an incorrect ex¬ 
pression for the precession of the longitude of periastron. I 
have found 20° as a lower limit for the inclination of the B 
star with respect to the orbital plane which does not change 
the conclusions concerning the neutron-star birth kick in this 
system given in Kaspi et al. (1996). 
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